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We present the diagrammatic technique for calculating the free energy of the matrix eigenvalue 

I I model (the model with arbitrary power (3 by the Vandermonde determinant) to all orders of \/N 

^ . expansion in the case where the limiting eigenvalue distribution spans arbitrary (but fixed) number 

^^ ' of disjoint intervals (curves). 

o. 

O ■ 1 Introduction 

o 



Exact solutions of various matrix models have been considerably put forward recently. The progress 

in finding asymptotic expansions is mostly due to a geometrization of the picture. Attained first for 

a mere large- A^ limit of the Hermitian one matrix model (IMM) |27j, ^1] ^J, 01 , it was almost 

simultaneously put forward for solutions of two-matrix Hermitian matrix model (2MM) ^21 ■, |1SJ , |2^ , 

[1], [3]. Combined with the technique of the loop equation [HI, which is the generating function for 

^ ' Virasoro conditions in matrix models U2]) it permitted to generalize the moment technique of [2] to 

I> ■ finding the subleading order correction to the IMM free energy first in the two-cut case ^ and then 

k> I in the multicut case [2^]) ESI) |Sj- Almost simultaneously, the same correction was found in the 2MM 

5—1 \ case ITO] and in the case of normal matrices (see (20] , [2B! ) • 

Next step in constructing asymptotic expansions pertains to introducing the diagrammatic tech- 
nique describing terms of expansion in the way very similar to the one in (quantum) field theories: the 
expansion order is related to the number of loops in the diagrams, and n-point correlation functions 
as well as the free energy itself are presented by a finite sum of diagrams in each given order of the 
expansion. This technique, first elaborated in [TBI for correlation functions in the IMM case, was then 
developed for pure and mixed correlation functions in the 2MM case |21) . |22) and then turned into 
a complete solution for the free energy first in the IMM case ^, and then, eventually, in the 2MM 
case 0. 
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We expect that generalizations of the diagrammatic technique are a very power tool for investi- 
gating various niatrix-model-like problems. We address the eigenvalue model, which describes a gas 
of N particles dwelling in the potential field and having anyonic statistics expressed by a power of 
the Vandermonde determinant in coordinates of the particles. This is the Dyson (or Laughlin) gas 
system, and we investigate a one-dimensional case of this system in this paper. (The first investigation 
in two-dimensional case was done recently in |^ where two first subleading order corrections were 
derived but without involving a diagrammatic technique.) We elaborate the diagrammatic technique 
for consistent calculation of corrections of all orders and find separately the only correction term that 
cannot be produced from this technique. 

The paper is organized as follows: in Sec. [21 we formulate the problem and present the loop 
equation together with necessary definitions. In Sec. |31 we introduce elements of our diagrammatic 
technique for calculating resolvents (loop means), and we summarize this technique in Sec. |1J In the 
next section El we invert the action of the loop insertion operator and obtain the free energy of the 
model, presenting first few corrections. This enables us to obtain all correction terms except two: 
the first one is just a subleading term in the IMM case (see jg]), while the second term has not been 
known previously. We calculate it separately in Sec. 

2 Eigenvalue models in l/A^^-expansion 

Our aim is to show how the technique of Feynman graph expansion elaborated in the case of Hermitian 
one-matrix model JH]) |Z] can be applied to solving the (formal) eigenvalue model with the action 



N 

n 



dx.A^^{x)e-t^'"''"'''=e-^, (1) 



where V{x) = X^n>o ^n^;" and 



h = ^ (2) 

is a formal expansion parameter. The integration in (^ goes over A^ variables Xj having the sense 
of eigenvalues of the Hermitian matrices for /3 = 1, orthogonal matrices for f3 = 1/2, and symplectic 
matrices for /3 = 2. In what follows, we set /3 to be arbitrary positive number. The integration may 
go over curves in the complex plane of each of N variables Xj. For (3 ^ 1, no topological expansion 
in even powers of h exists and we rather have the expansion in all integer powers of h. Customarily, 
to = hN is the scaled number of eigenvalues. We also assume the potential V{p) to be a polynomial 
of the fixed degree m -|- 1, or, in a more general setting, it suffices to demand the derivative V'{p) to 
be a rational function |2()j . 

The averages corresponding to partition function (^ are defined in a standard way: 



(Qix)) 



i / DX Q{X) exp f -iV^tr V{X)] , (3) 



and we introduce their formal generating functionals: the one-point resolvent 



mp) = nvr^t^ (4) 



fc=0 



as well as the s-point resolvents {s > 2) 



»F(p,....p,) = ^(-^j Y. 



A:i,...,fc,=l ^1 ^^ 



2-s 



( hY~' I 1 1 



^vjJ V-VT^''''^^^^ <^' 



conn 



where the subscript "conn" pertains to the connected part. 

These resolvents are obtained from the free energy T through the action 



W(j)\....,Ps) = -^ 



2 d d dr 



d d d 

W^(pi), (6) 



dv{ps) dv{p,_i) dv{p2 

of the loop insertion operator 

^- = -V^— (7) 

Therefore, if one knows exactly the one-point resolvent for arbitrary potential, all multi-point resolvents 
can be calculated by induction. In the above normalization, the /i-expansion has the form 

oo 

W{pi,...,Ps)=Y,h'W,/2{pi,...,Ps). S>1, (8) 

r=0 

where it is customarily assumed that it corresponds in a vague sense to the genus expansion in the 
usual Hermitian models with possible half-integer contributions. It is often written as a sum over 
all — integer and half-integer — g = r/2. 

The first in the chain of the loop equations of the eigenvalue model (^ is 

/cd 2mp-uj V / 

Here and hereafter, Cx> is a contour encircling clockwise all singular points (cuts) of W{lo), but not 
the point to = p; this contour integration acts as the projection operator extracting negative part of 
V'{p)W{p). Using Eq. Q, one can express the third second term in the r.h.s. of loop equation ^ 
through W{p), and Eq. Q becomes an equation on one-point resolvent Q. 

The /3-dependence enters @ only through the combination 



7 = V^-V^, (10) 

and, assuming /3 ~ 0(1), we have the free energy expansion of the form 

oo oo 

:f ^ T{h, 7, to, ti, i2, . . .) = E E ^''+'' V-^fe.i- (11) 

fc=0 1=0 



Substituting expansion (jSJ in Eq. ©, we find that Wg{p) for g > 1/2 satisfy the equation 

9-1/2 

K - 2W0{P)) Wg{p) = Y. Wg>{p)W,^g,{p) + ^^^Wg^,{p)+J-Wg^y2{p), ^ 

where K is the hnear integral operator 

In Eq. (|12p . Wg(p) is expressed through only the Wg.{p) for which gi < g. This fact permits developing 
the iterative procedure. 

In analogy with Hll|) . it is convenient to expand multiresolvents Wg{-) in 7: 

[9] 

Wgipi, ...,p,) = Y^ J^a-2k Wk,2g-2k{PU ■■■,Ps) (14) 

fc=0 

Then, obviously, ()12() becomes 

^ - 21^o,o(p)) WfcXp) = E Wk,Mip)Wk-k,,i-iAp) + ^^TTT^'"'''^^^ + 7^W^fe,z-i(p). (15) 

fcl>0,ii>0 ^"' " 

fcl+!l>0 

The form of loop equation Q is based exclusively on the reparameterization invariance of the 
matrix integral, which retains independently on the details of eigenvalue density distribution. We 
assume that as A^ — > 00, the eigenvalues fill in some segments in complex plane, dependently on the 
shape of potential V{X). For polynomial potentials, the number of segments is finite and the contour 
C-jy of integration in (|12j) encircles a finite number n of disjoint intervals 

n 

^= |J['"2i-l>/"2i]> /^l < /^2 < • • • < /"2n- (16) 

i=l 

Recall that 



all Wk^i{p) are total derivatives, 



Wk,i{p)\p^oc = -SkAfi + 0{l/p^), (17) 

P 



Wk,i{p) = -g^f^u k,l>0. (18) 



The solution in the large- A^ limit coincides with the solution of the Hermitian one-matrix model and 
satisfies the equation (in fact, the above normalization was chosen in a way to ensure the coincidence 
of this equation with the one in the IMM case) 

^^^'^^^Wo{X) = {Wo(p)f. (19) 



c„ 2TTip- X 



Actually, Wo{p) must be Wo,o(p) in our notation; we however preserve the old notation assuming this 
identification in what follows (for shortening the presentation). 

Recall the solution of Eq. (|19() . Deforming the contour in Eq. (|19|) to infinity, we obtain 

{Woip)f = V'ip)Woip)+i ^^^Wo{X), (20) 

where the last term in the r.h.s. is a polynomial Pm-i of degree m — 1, and the solution to H2U() is 
then 

Woip) = ^Vip) - lv^'(p)2 + 4P^_,(p) ^ ly'(p) _ y(p), (21) 

where the minus sign is chosen in order to fulfill asymptotic condition p7j) whereas the function y{p) 
is defined as follows. For a polynomial potential V of degree m + 1, the resolvent VFo(p) is a function 
on complex plane with n < m cuts, or on a hyperelliptic curve y^ = V'{p)'^ + 4Pm-i(p) of genus 
g = n — 1. For generic potential V^(-'l^) with m — > oo, this curve may have an infinite degree, but we 
can still consider solutions with a finite genus, where a fixed number n of cuts are filled by eigenvalues. 
For this, we separate the smooth part of the curve introducing 

y{p)^M{p)y{p), and f{p) = \\ (p - fia) (22) 

with all branching points /x^ distinct. The variable y defines therefore the new, reduced Riemann 
surface, which plays a fundamental role in our construction. In what follows, we still assume M{p) to 
be a polynomial of degree m — n, keeping in mind that n is always finite and fixed, while m > n can 
be chosen arbitrarily large. 

From now on, we distinguish between images of the infinity at two sheets — physical and unphysical — 
of hyperelliptic Riemann surface (|22jl respectively denoting them oo+ and oo_. We often distinguish 
between variables of the physical and unphysical sheets placing the bar over the latter. By convention, 

we set y|p_^oo+ ~ p", and M{p) is then-"^ 

1 V'(w) 

M{p) = --resoo+du;- ^-^. (23) 

2 [w — p)y{w) 

Inserting this solution in Eq. (|21l) and deforming the contour back, we obtain the planar one-point 
resolvent with an n-cut structure, 

2 Jct, 27rz p - a y(\) 
Let us now discuss the parameter counting. We introduce the filling fractions 

where Ai, i = 1, ... ,n — 1 is the basis of ^-cycles on the hyperelliptic Riemann surface (|22j) (we 
may conveniently choose them to be the first n — 1 cuts). Adding the (normalized) total number of 
eigenvalues 

f y(A)— = reSoo+y(A)(iA (26) 



to 

'Ct, 



^By a standard convention, resaadx/x = — 1, and the direction of the integration contour at the infinity point therefore 
coincides with the direction of contour for integrals over C-d and over the set of yl-cycles, see below. 



to the set of Si, we obtain n parameters, to which we add, foUowing ((T7|l. the asymptotic conditions 

1 / dX \''V'(X) , 
-'"'"-- 2 t 2^,^- * = ».■■■'"■ (27) 

In this article we consider fiUing fractions (|25|) as independent parameters of the theory. We need 
this assumption to interpret random matrix integrals as generating functions of discrete surfaces. 
Other assumptions are possible, but we do not consider them in this paper. In other words, we 
consider only the perturbative part of the matrix integral, and the filling fractions are fixed because 
the jumps between different cuts are non-perturbative corrections in h. In particular, this imposes 
restrictions 

^5.=0, , = 1,...,„-1, g4yto = 0. (28) 

That implies that, for k + I > 0: 

/ WkM^P2,...,Ps)=0 (29) 

In addition, we impose another assumption. The zeroes bj of M{p) are called double points. In 
some sense, they can be considered as cuts of vanishing size. We require that those degenerate cuts 
contain no eigenvalue to any order in the h expansion. We therefore demand 

/ Wk^i{i,p2,...,Ps) = ^ (30) 

for any contour which encircles a zero of M. 

This assumption, together with loop equation (jlSj) . suffices for proving that for every /c, I, s, except 
{k, I, s) = (0, 0, 1) and {k, I, s) = (0, 0, 2), the function Wk^iipi, ■ ■ ■ ,Ps) has singularities in the physical 
sheet only at the branch points //„. In particular, it has no singularities at the double points in the 
physical sheet. 

3 Calculating resolvents. Diagrammatic technique 

In this section, we derive the diagrammatic technique for model (^, which is a generalization of 
technique in ^S] , [7j . Our main goal is to invert loop equation (fT^ to obtain the expression for Wk (p) 
for any k > 1/2. 

3.1 A piece of Riemann geometry 

The main notion is again the Bergmann kernel,'^ which is the unique bi-differential on a Riemann 
surface T,g that is symmetrical in its arguments P,Q G T,g and has the only singularity (a double pole) 
at the coinciding arguments where, in any local coordinate r, it has the behavior (see |21], j23j ) 



^It is a double derivative of the logarithm of the Prime form. 



with Sb{P) the Bergmann projective connection associated to the local coordinate r. We fix the 
normalization claiming vanishing all the integrals over yl-cycles of B{P, Q): 

i B{P, Q) =0, fovi = l,...,g. (32) 

We then have the standard Ranch variational formulas relating B(P, Q) with other objects on a 
(general, not necessarily hyperelliptic) Riemann surface: 

^ B{P, Q) = 2B{P, [^«])5(K], Q), (33) 



and 

(p B{P,Q) = 2mdwi{P), (34) 

where //q, is any simple branching point of the complex structure. Then, by definition, in the vicinity 
of /ia, 

BiP, g)|Q-.^, = B{P, [fi^]) (^7=^ + 0{./^^^)dq^ , (35) 

and dwi{P) are canonically normalized holomorphic differentials: 

dwiiP) = 5ij. (36) 

Besides these formulas, we need another, rather obvious, relation: for any meromorphic function 
/ on the curve, we have: 

d/(P) = ^jf B{P,OfiO (37) 

where the contour Cp encircles the point P only, and not the poles of /. 

We also introduce the 1-form dEg^g^^P), which is the primitive of B{P,Q): 

dEQ,,,{P)=j B{P,C), dEQ,,,{P)\p^Q= /^^(-^) + finite. (38) 

Then, obviously. 



/ dEQ,g,{P)=0. (39) 



The form dEQ^q^{P) is single- valued w.r.t. P on the Riemann surface and multiple- valued w.r.t. the 
variable Q: from (|34j) . 

di?Q+x ,,„(P) = 2T,idw,{P)+dEQ,,,{P), 



^B, 



dEQ+,,AP)=dEQ^,,{P). 



^A, 



where the reference point go will happen to play no role and it must disappear eventually from all 
formulae. Technically it is convenient to have go ™- the non-physical sheet. 



We can now express the 2-point resolvent Wo{p, q) in terms of B{P, Q). We let p and p denote the 
complex coordinates of points on the respective physical and unphysical sheets. Then, 

dV'ip) dpdq 

dpdq——— = -B{p,q)-B{p,q) - 



dV{q) ' ' {p — qY 

since it has double poles with unit quadratic residues &t p = q and p = q. The 2-point resolvent (|21j) 
is nonsingular at coinciding points; therefore, 

dpdq^^= -l{B{p,q) - B{p,q)\ (40) 

and 

3.2 Inverting the operator K — 2Wq{p) 

We can determine corrections in h iteratively by inverting loop equation H12|) . All multi-point resolvents 
of the same order can be obtained from Wg{p) merely by applying the loop insertion operator ^y^ ■ , . 

In the IMM case, a natural restriction imposed on the free energy is that all the higher free energy 

(k) 
terms Fg must depend only on /Zq and a finite number of the moments Ma , which are derivatives of 

{k — l)th order of the polynomial M{p) at branching points, allowed no freedom of adding the terms 

depending only on tg and Si to J-g. In model (^, such a restriction cannot be literally imposed, as we 

demonstrate below, and we instead claim that 

di, (42) 



which was a consequence of locality in IMM and can be taken as a defining relation in the /3-model. 

The first step is to find the inverse of the operator K — 2Wo{p). It was found in I^, that if 
f{p) is a function whose only singularities in the physical sheet are cuts along D, which vanishes at 
oo like 0(l/p^) in the physical sheet and has vanishing A-cycle integrals, then, having dEq^q{p) = 

dEg,qoip) -dEg^q^ip), 

d£p,g{K - 2Wo{q))f{q) ^ ^. i ^E±M^(^K - 2Wo{q)).fiq) = fip) dp, (43) 

2tti Jc^ 2y{q) 

where y{q) = — 2 — W^o(?) is the function introduced in H21|) and integration contour lies in the physical 
sheet. Indeed, in the standard notation for the projections [(5((?)]± that segregate the respective 
polynomial and regular parts of Q{q) as g — > 00, we have {K — 2Wo{q)).f{q) = 2[y(g)/(g)]_ = 
2y{q)f{q) — P{q) where P{q) is a polynomial [y(g)/((?)]+ of degree degF' — 2, and thus: 

Let us compute the part of the integral involving the polynomial P{q)'- 



/ dEg^^ip) dq ^ + 2m V Res dEq^q{p) dq ^^ 
/ dEg^g{p)dq^ 

Jc^ yw 



= 0. (45) 

The first equality is obtained by changing the name of the variable q ^>- q with accounting for P{q) = 
P{q), y{q) = —y{q)- The second equality comes from deforming the contour Cx> to —C-p picking 
residues at branch points. The third equality holds because P{q) is a polynomial and thus has no 
singularity at branch points, whereas zeros of y{q) are canceled by those of dEq^q{p), so the residues 
vanish. 

Therefore we have: 

= ^ / , {dEg,goip) + dw{p)) dq f{q) 

= ^ ^, dEq^q^{p)dqf{q) 

= - Res dEq^qo{p)dqf{q) 

= f{p)dp. (46) 

The second equality holds because dEq^qf^{p) has no singularity at q ^^ p and we can push the inte- 
gration contour for q to infinity (in the physical sheet), which gives zero. In the third equality, the 
contour C'jy is a contour which encloses the ^-cycles and V. When we cross the cycle Ai, dEq^q^^ijp) 
jumps by the corresponding holomorphic differential dwi{p). The fourth equality holds because dwi{p) 
is independent of g, whereas the integral of f{q) along any yl— cycle vanishes by our assumption (|29|) 
(this is nothing but the Riemann bilinear identity). Then the contour C^ is deformed in the physical 
sheet, into a contour which encloses only the point p (recall that we have assumed that f{q) has no 
singularity in the physical sheet, and vanishes like 0{l/q^) at cxd), and the final result comes from the 
fact that dE has a simple pole (see eq. (|5H)) ). 

This proves that for all functions of the type W^^u the integration i^ §^ — ''iil '' ^^^^ ^^^^ ^^^ 
inverse of if — 2Wo{q). 

Notice that if f{q) is a function which has poles only at branch points fia or at double points bj, 
by moving the integration contours we have: 

^ / iEiM^ /„) . Y, Res 'A^ m + E R? ^%rr^ (/(«) - /(?»■ («) 

^7r« Jc-r, ^y{Q) -fi^'"'^ 2y(g) , i^^j 2y(9) 



2m 



This relation provided a basis for the diagrammatic representation for resolvents in IMM ^^, and 
we show that it works in the case of /3-model as well. Let us represent the form dEq^q[p) as the vector 
directed from p to q^ the three-point vertex as the dot at which we assume the integration over g, 

dq 1 

27ri 2y(q) ' 



§C 2m 2 ( ) ' ^^^ ^^^ Bergmann 2-form B{p,q) as a nonarrowed edge connecting points p and q. 



Let us also introduce a new propagator dpdy{q) denoted by the dashed line. 

The graphic representation for a solution of (|12|) then looks as follows. We represent the multire- 



solvent Wg'{pi,... ,ps 
derivative -^Wg'ipi 



as the block with s external legs and with the index g' . We also present the 
■ ■ ,Ps) as the block with s + 1 external legs, one of which is the dashed leg that 




starts at the same vertex as pi. That is, we obtain (cf. ^^) 



HZ) 



(48) 

which provides the diagrammatic representation for Wk{pi, ■ ■ ■ ,Ps)- 

Recall the diagrammatic formulation of IMM (7 = 0). There the multiresolvent Wkfi{pi, . . . ,ps) 
can be presented as a finite sum of all possible connected graphs with k loops and s external legs and 
with only three- valent internal vertices (the total number of edges is then 2s + 3A; — 3, and we assume 
s > 1 for A; > 1 and s > 3 for k = 0) and such that in each graph we segregate a maximal rooted tree 
subgraph with all arrows directed from the root. This subtree comprises exactly 2k + s — 2 arrowed 
edges. We then choose one of the external legs, say, pi (the choice is arbitrary due to the symmetry 
of Wkfiipi, ... ,Ps)), to he the root vertex the tree starts with; for each three- valent vertex there must 
exist exactly one incoming edge of the tree subgraph. All external edges (except the root edge) are lines 
corresponding to B{p,q) and are therefore nonarrowed. This subtree therefore establishes a partial 
ordering of vertices: we say that vertex A precedes vertex B if there exists a directed path in the 
subtree from A to B. Internal nonarrowed edges are again B{r, q) but we allow only those nonarrowed 
edges for which the endpoints, r and q, are comparable. If r = g, then, for the tadpole subgraph, we 
set B{r,f), where f is the point on the other, nonphysical sheet. At each internal vertex, denoted by 
• , we have the integration f,(g) 2^ JTT' ^^i^^ the arrangement of the integration contours at different 
internal vertices is prescribed by the arrowed subtree: the closer is a vertex to the root, the more outer 
is the integration contour.^ 



3.3 Acting by spatial derivative 

As a warm-up example, let us consider the action of the spatial derivative d/dpi on Wkfi{pi, • • • ,Ps)- 
We place the starting point of the ( fictitious) dashed directed edge to the same point pi and associate 
just dx with this starting point. Recall that the first object (on which the derivative actually acts) 
is dE^^fj{pi), then comes the vertex with the integration ^ f^M ^jt^y; then the rest of the diagram, 

which we denote as F(ri). We can present the action of the derivative via the contour integral around 
pi with the kernel B{pi,S,)'. 



dpi 



d I / dEr^^f,ipi)dr] 
dpi l^J^t,"' 27rz y{rj) dpi 



Fill) 



Res 

?^pi Jc. 



B{pi,OdEr,,fj{0 dr] 



(v) 



d^ 



2m y{ri) 



Fill), 



(49) 



Since no propagator connects noncomparable vertices, the mutual ordering of the corresponding integration contours 
is irrelevant. 
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where pi lies outside the integration contour for rj. The integral over ^ is nonsingular at infinity, so 
we can deform the integration contour from Cp^ to C^ > Cj^ . 



2-Ki d^ JqM 2-Ki y{if) Tc''^^>c^'^^ ^vri d^ 27ri y{rj) 

We now push the contour for ^ through the contour for rj, picking residues at the poles in (^, at ^ = r], 
^ = fj and at the branch points. We then obtain 

YRe. ^'"'-^'f-"'^' "", , F(„) - /, , |^F(,) + /, _ I^F(„) (51) 

where the first integral experiences only simple poles at the branching points. The residue is unchanged 
by using I'Hopital rule, replacing B{^,pi)/dy{^) by dEcc{pi)/2y{S,), that is, we have 

\-j>^,, dE^,M)dy{OdE,,,{0 dr^ r Bjpuv) - Bjpuv) 

C'^'^T^^^- '^y^^di 2TTiy{r]) ^ '^ Jc^ 2^i y(r/) ^ '^ 

Pushing contour of integration for ^ around the branch points back through the contour for rj, we pick 
residues at .^ = t? and ^ = fj, which both give the same contribution, and we obtain 

dE^^^{pi)dy{OdEr,,fj{0 dr] ^ 
d«>d'" 2TTi2y{Od^ 27riy(r?) ^'^^ 



T> '^'-D 



^ dE,,,{p,)dy{v) dy _r B{p,,r,)-B{p^,fi) 

Jc!^^ y{v) drj 27ri y(ri) Jc^ 27ri y{r]) 

We evaluate the last two terms by parts, so it remains only 

dE^^^ipi)dy{OdE^,fj{() dv p^.^I dE,,f^{p,)dri OF JT]) 



C«'>4'') 2TTi2y{Od( 2TTiy{r]) JcM 27Ti y{i]) di] 



We have thus found that 



dpi y/c*,"' 27r« y{r]) dpi J 

dE^^^{pi)dy{OdE^,^{i) dy p( . I dE^^^jpi) drj dFjr,) 
C«)>c(,'') 27ri2y{OdC 2TTi y{7]) ^'^> fc^ 2t: i y{r,) Oy ' 



and we can graphically present the action of the derivative as follows: 

dp,^[F{r,)] = Pi. ^^ ' ^ •{f{v)}+ Pi^{|F(r?)} 



(55) 



*Here, as in 0, comparison of contours is equivalent to their inside/outside ordering. 
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We therefore see that using relation H55() we can push the differentiation along the arrowed edges 
of a graph. It remains to determine the action of the derivative on internal nonarrowed edges. But for 
these edges (since it has two ends), having the term with derivative from the one side, we necessarily 
come also to the term with the derivative from the other side; combining these terms, we obtain 



dpB{p,q)+dgB{p,q) = f 

JC-n 



2TTi di 



and we can deform this contour to the sum of contours only around the branching points (sum of 
residues). Then we can again introduce y'{C)di/y{Cj and integrate out one of the Bergmann kernels 
(the one that is adjacent to the point q \i p > q ov p \i q > p; recall that, by condition, the points p 
and q must be comparable). 

That is, we have 



{dp + dq) P- 'Q = P' *i ■ g = P' 



(56) 

It becomes clear from the above that we must improve the diagrammatic technique of /3-model 
in comparison with the IMM by including the dashed lines (we do not call them propagators as 
they have a rather fictitious meaning); being treated as propagators, they however ensure the proper 
combinatorics of diagrams. Indeed, from ()55() and H57|) it follows that the derivative action on the 
"beginning" of the dashed line is null, dpdp = 0, and when this derivative acts on the "end" of this 
line, we merely have 

dW{Q)=y"{Q), (57) 

which we denote symbolically as two dashed propagators ending at the same vertex. If we continue to 
act by derivatives d/dq, then, obviously, when k dashed propagators are terminated at the same vertex, 
we have the kth order derivative y^ ' (q) corresponding to them. We then already have three variants 
of incorporating dashed lines into the play. Recall that the dashed lines, as the solid nonarrowed 
lines, can connect only comparable vertices (maybe, the same vertex) and the starting vertex must 
necessarily precede the terminating vertex (they may coincide). 

The first case is where we have two adjacent solid lines and exactly one outgoing dashed line. Then, 
no incoming dashed lines or extra outgoing dashed lines is possible and the both solid lines must be 
arrowed: one is pointed inward and the other outward. We also distinguish this case by labeling the 
corresponding vertex by the white dot: 



dq 
cM 2m y{q) ' 

(58) 

The second case is where we have two adjacent solid lines (one of them is necessarily incoming 
directed line, but the other can be either directed outward or be nonarrowed internal or external line), 
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no outgoing dashed line, and k > 1 incoming dashed hnes (this means that we must have at least k 
white vertices preceding this vertex in the total graph). In this case, we denote the corresponding 
vertex by a solid dot: 

k 

q ^ ^ JcM 27ri y{q) 

(59) 

The last case is when we have just one (incoming arrowed) solid line adjacent to the vertex. We 
then have exactly one outgoing dashed line and k + 1 incoming dashed lines (A; > 0). In this case, we 
also denote the vertex by the solid dot: 

k 

y^^+^\q)dq 

(60) 

External lines of Wk,i{pi, . . . ,ps) are either the root vertex pi or nonarrowed propagators B(q,pi); 
no external dashed lines are possible. 

Our general rule for assigning white and black colors to vertices is as follows: if there is no factors 
y^ '{q) standing by the vertex, it is white; if there are such factors, the vertex is painted black. 

Now comes the question what happens if we act by the loop insertion operator d/dV{p) on the 
elements of the constructed diagrammatic technique. 

3.4 Acting by the loop insertion operator 

We now extend our diagrammatic technique by incorporating the action of loop insertion operator 
((7|) on its elements. The action of the loop insertion operator on the Bergmann differential and its 
primitive was presented in [7j, so we do not describe it here. We can graphically present the action of 
d/dV{r) as 



'-^'' ~ i«^^' *^°' 



j^p 



aV{r) ^' -^^- ^'--^^r-.^^-, gy(^) 



-k!)-^P, ?5W7) ^' 'P = ^-*0-^P = 9— CH- 



e e e 



(61) 



where in the first case we must also take into account the variation of the y{p) factor in the denominator 
of the measure of integration in p at the right vertex (irrespectively, white or black), and the proper 
contour ordering is assumed. All the appearing vertices are of white color as they do not contain 
additional factors of type y^^'H)- In the second case, it is our choice on which of edges to set the 
arrow. Recall, however, that the points P and Q were already ordered, as prescribed by the diagram 
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technique. That is, if "P > Q", we must choose the first variant and if "Q > P", we must choose the 
second variant of arrows arrangement in order to preserve this prescription. 

We now calculate the action of d/dV{r) on the dashed propagator. Obviously, 

but any attempt to simplify this expression or to reduce it to a combination of previously introduced 
diagrammatic elements fails. This means that we must consider it a new element of the diagrammatic 
technique. With a little abuse of notation, we visualize it by preserving k dashed arrows still landed 
at the vertex "(7" with the added nonarrowed solid line (third or second, depending on whether this 
vertex was of the second ()59|) or third H60() kind) corresponding to the propagator B{r,q). The vertex 
then change the coloring from black to white because it contains y^ '{q) factors no more. We then 
represent graphically these vertices as 



av(r) 



k 



\ / 




P^^^i^) = P^^^i^)^ U2^)W''^'^'^^ '-'■ 



V 



dq d^ 

^n>{r,q), k ^ I. 

(63) 



and 



av{r) 




f da 9^^+^ 

(64) 
for the respective cases (|^ and (|Hn|) . Recall that here necessarily r > q. 

The reason why we prefer to keep arrowed propagators entering "white" vertices will be clear after 
the next step when we consider the subsequent application of the spatial derivative d/dq to the new 
object -§^B{r,q). 

Recall that the very application of this operator assumes that we have somewhere in the preceding 
vertices a vertex of type H58() . i.e., the vertex with a new white vertex. Then, successively applying 
derivative when moving up from the root over branches of the tree subgraph, we have two possibilities. 
The first case is where the vertex "r" is an external vertex of the graph. Then, when we reach the 
vertex "g" we just set the extra derivative 

on the corresponding Bergmann kernel, which corresponds to adding an extra incoming dashed line to 
either the case (|H5]) or to the case (|64p . depending on the type of the "white" vertex. The second and 
most involved case occurs when the vertex "r" is an internal vertex of the graph. Then, the derivative 
must successively act on both its ends, and we want to express the action on the end "r" in terms of 
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the diagram technique constructed above. We have from 1)56^ 



(65) 

and it remains to act by k derivatives d/dq using the rules formulated in (|55|) . We then produce a 
number of diagrams, but all of them will be of one of the type indicated above. It becomes clear why 
we prefer to preserve incoming dashed lines in notation (jHSJ and (|(i4)) . When applying derivatives d/dq 
in the r.h.s. of H65|) a part of these dashed lines that does not act on the last nonarrowed propagator 
B{r,^) appear again as the derivatives y^^'{xj) at intermediate vertices "x/' (^ < Xj < 9 < ''')• 

The very last step in constructing the diagram technique is to consider the action of the loop 
insertion operator on -0-^B{q,r). Then, applying now relation (|6H) . we obtain 



dV{p)dq^^''''" ''^^^S,'^ 



a d^ 

:B{r, 



(66) 



and we observe the appearance of the last remaining structure of the diagram technique — the ( "white" ) 
vertex "^" at which a number s < A; of dashed lines terminate and which have two incident nonarrowed 
solid lines corresponding to the propagators B[^,p) and B{^,r). This vertex generalizes vertex H63|) 
in a sense that the action of derivatives must be now distributed among these two nonarrowed lines 
(see (jMI below). 

It may appear the situation where we have the closed solid loop (the propagator B{^,^) as shown 
below) 







(67) 

In this case, however, as soon as A; > 0, the derivatives act on both ends of the propagator S(^, ^), 
and we can use (|65() to distribute them. Therefore, such a diagram (with closed loop of S-propagator) 
enters the diagram technique only in the case k = 0. 

We are now ready to present the complete diagram technique for multiresolvents of the /3-model. 
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4 Feynman diagram rules 

We therefore have the following components of the diagrammatic technique. We absorb all the factors 
related to Bergmann kernels B{p,q) and dEq^q{p) and all the factors with y{^) and its derivatives 
into vertices in accordance with the rule: we associate to a vertex the solid arrowed propagator 
that terminates at this vertex (recall there is exactly one such propagator) and associate nonarrowed 
propagator B{p, q) or its derivatives with the vertex that correspond to the minimal variable among p 
and q. The ordering of vertices is implied from left to right (as will be assumed in most of appearances 
below). 

The vertices with three adjacent solid lines 



k 



's¥i ~ i.«''^«.^-<''>d!(i)|^(''('''«^<''-«)- 's»^ 



p^/i 



r 



(68) 



Here, by construction, ^ < r and ^ < p, or r and/or p can be external vertices. If q here is an external 
vertex, then A; = and r and p must be also external vertices; 



-g? - li/''^^^^2^f^^^~^^^ ^^0. 



(69) 



The vertex q can be external. 

k 



q^^p ^ l,,dE,,,[q),^--^—^Bir,0. k > 0, r > C. 



d^ d'' 

(70) 
Here r can be external vertex. If q is an external vertex, then A; = and r is also an external vertex. 



/ d^ 



(71) 
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Here q can be an external vertex (and r is not by condition). 



^^f^P ~ ^(O'^^?'?^^) 



d( 



2m y{^)' 



(72) 



Here q can be an external vertex. 

The vertices with two adjacent solid lines 

k 

pSy' ~ I dE,,,{pf^^^^B{r,q), k>l. 

yq /c<,«' 27ri y{q) 

(73) 

Here by construction q < r and r can be an external vertex, while p cannot be an external vertex. 

k 

q r J^M 2m y{q) 

(74) 

Here p cannot be an external vertex. 

k 

dq d^+^ 



s /■ 



'^/ ^ ^/c(?)'''-"^^)^^^5?TT^(-'^)' ^>0. 



(75) 



Here always (7 < r or r is an external vertex. If p is an external vertex, then A; = and r is also an 
external vertex. 



P—kA^t ~ J (p dEq^qip)- 



Tc!^) "'''^'2my{q) 

(76) 
Here p can be an external vertex. 
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p^^g ~ 7 t,, dE,,ip)^-^, k > 0. 



The vertex with one adjacent solid line 

k 

(77) 

If p here is external, then k = 0. 

When calculating Wg{pi, . . . ,ps) we take the sum over all possible graphs with one external leg 
dEq^q{pi) and all other external legs to be B{^j,pi) and, possibly, their derivatives w.r.t. internal 
variables ^j in accordance with the above rules such that arrowed propagators constitute a maximum 
directed tree subgraph with the root at pi and the bold nonarrowed lines can connect only comparable 
vertices (and may have derivatives only at their inner endpoints) and arrowed dashed lines can also 
connect only comparable vertices (the arrow is then directed along the arrows on the tree subgraph). 
The diagrams enter the sum with the standard symmetry coefficients. 

In full analogy with the IMM case, the order of integration contours is prescribed by the order of 
vertices in the subtree: the closer is the vertex to the root, the more outer is the integration contour. 
In contrast to the IMM case, the integration cannot be reduced to taking residues at the branching 
points only; all internal integrations can be nevertheless reduced to sums of residues, but these sums 
may now include residues at zeros of the additional polynomial M{p) on the nonphysical sheet and, 
possibly, at the point oo_. 

In complete analogy with the IMM case, in the next section, we use the iif-operator introduced 
in j2j in order to invert the action of the loop insertion operator and obtain the expression for the free 
energy itself. 

5 Inverting the loop insertion operator. Free energy 

5.1 The if -operator 

We now use the operator that is in a sense inverse to loop insertion operator ((J)) and was introduced 
in [7] in the IMM case. It has the form^ 

H- = -Tesoo+V{x) ■ --res^_V{x) ■ -to • - ^5, ^ •. (78) 

2 = 1 ^ 

The arrangement of the integration contours see in Fig. 1. We now calculate the action of H on 
the Bergmann bidifferential B{x,q) using again the Riemann bilinear identities. We first note that 
B{x,q) = {dxdEx^qQ{q)) dx and we can evaluate residues at infinities by parts. Then, since dEx^qgiq) is 
regular at infinities, for V'{x) we substitute 2y[x) + 2io/x as x — > oo-|_ and —2y{x) + 2to/x as x — > cxd_ 
thus obtaining 

-reSoo+ ( y{x) + — j dEx,qo{q)dx + reSoo_ ( -y{x) + — j dEx,qo{q)dx 



^This definition works well when acting on l-forms regular at infinities. Otherwise (say, in the case of Wo(j>)), the 
integral in the third term must be regularized, e.g., by replacing it by the contour integral around the logarithmic cut 
stretched between two infinities. 
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Figure 1: The arrangement of integration contours on the Riemann surface. 



n-l 

1=00+ 



=1 JB. 



todE^,g,{q) _ -Y^Sii B{q,x), (79) 



whence the cancelation of terms containing to is obvious, and it remains to take the combination 
of residues at infinities involving y{x). For this, we cut the surface along A- and S-cycles taking 
into account the residue at x = q. The boundary integrals on two sides of the cut at Bi then 
differ by dEx^qoiq) — dE^_^_i q(^{q) = 0, while the integrals on two sides of the cut at Ai differ by 

dEx^qg{q) — dE^^i ^q^iq) = §^, B{q,x), and we obtain for the boundary term the expression 



n— 1 „ „ 

Y^ i y{x)dx j> B{q,C), 



which exactly cancel the last term in ()79() . It remains only the contribution from the pole at x = g, 
which is just —y{q). We have therefore proved that 

H-B{;q) = -y{q)dq. (80) 

Let us now consider the action of H on Wk,i{-) subsequently evaluating the action of loop insertion 
operator (jT)) on the result. Note first that the only result of action of d/dV(p) on the operator H 
itself are derivatives dV{x)/dV{p) = — l/(p — x) (and recall that by definition \p\ > \x\, i.e., instead 
of evaluating residues at infinities one should take residues at x = p, and we obtain 

^^ {H ■ Wk,i{-)) = Wk,i{p) + H ■ Wk,i{;P)- (81) 

For the second term, due to the symmetry of Wk^i{p,q), we may choose the point p to be the root of 
the tree subgraphs. Then, the operator H always acts on B{-,^) (or, possibly, on its derivatives w.r.t. 
^) where ^ are integration variables of internal vertices. 

Let us recall the action of d/dV{q) on the elements of the Feynman diagram technique in Sec. EI 
Here we have three different cases. 

• When acting on the arrowed propagator followed by a (white or black) vertex, we use the first 
relation in ()61|). 



When acting on nonarrowed internal propagator ■4-^B{p, q), p > q, k > 0, we apply relation (1661) 



ak 



cik 

without subsequent representing the action of the derivative #x ^^ & sum of diagrams. We have 
no external i?-lines as we act on the one-loop resolvent. 
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• Eventually, when acting on dashed lines coming to a black vertex, using relation (|62() we obtain 
expression in H63|) : the action on dashed lines coming to a white vertex is null. 

We now consider the inverse action of the //-operator in all three cases. 

In the first case where it exists an outgoing arrowed propagator dEp^qg{^) (we can have only one 
such arrowed propagator as one line is external), then we can push the integration contour for ^ 
through the one for p; the only contribution conies from the pole at ^ = p (with the opposite sign due 
to the choice of contour directions in Fig. 1. We then obtain the following graphical representation 
for the action of the operator H in the first case: 

H- 



•—6-*% 



P = - Q' — >%P 
^ (82) 

In the second case, the vertex ^ in ()66() is an innermost vertex (i.e., there is no arrowed edges 
coming out of it). The 1-form y{(,)dS^ arising under the action of H (|8fl|l cancels the corresponding 
form in the integration expression, and the residue vanishes being nonsingular at the branching point. 
Graphically, we have 



H- 



^ q ^1 



0, k>0. 



(83) 



Eventually, in the third case, the inversion is rather easy to produce. Indeed, action of //-operator 
just erases the new /3-propagator arising in H63|) simultaneously changing back the color of the vertex 
from white to black: 

k k 



V ^ N ^ 



p— •^(^) = p- »V (^) k>l. 

I//- 






(84) 



For Hq ■ Wk{q,p) = Hq ■ gyr) Wk (p) , we obtain that for each arrowed edge, on which the action of 
produces the new (white) vertex, the inverse action of Hq- gives the factor —1, on each nonarrowed 
edge, on which the action of (O produces the new vertex accordingly to H66|). the inverse action of Hq- 
just gives zero, and at each black vertex, at which the action of changes the color to white and 
adds a new /^-propagator, the inverse action of Hq- gives the factor +1. 

As the total number of arrowed edges coincides with the total number of vertices and the contri- 
butions of black vertices are opposite to the contributions of arrowed edges, the total factor on which 
the diagram is multiplied is exactly minus the number of white vertices, which is 2A; + / — 1 for any 
graph contributing to Wk^iip)- We then have 

Hq - Wk,iiq,p) = -C^k + / - i)Wk{p) 
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and, combining with l\Hl\i . we just obtain 



{H, ■ Wk,iiq)) = -(2 - 2A: - Ot^tt^^M' (85) 



and, since all the dependence on filling fractions and to is fixed by condition (|42|) . we conclude that 

This is our final answer for the free energy. It permits us to calculate all J-k,i except the contribution 
at A; = 1,^ = (torus approximation in the IMM) and the second-order correction in 7 (the term 
.7^0,2) ■ The term J^i^ was calculated by a direct integration in |S]. We devote the special section 
below for the calculation of the term .7-0,2 • AH other orders can be consistently calculated. For this, 

■> J^ y(s)ds 

we only introduce the new vertex at which we place the (nonlocal) integral term ^^(^) ^ ^ /c\ — • 
In the IMM case, although this term was also nonlocal, it was possible to shift the starting point at 
the branching point Ha in the vicinity of every branching point; here it is no more the case and we 
must consider global integrations. Note, however, that it is only for the very last integration for which 
we must introduce nonlocal terms; all internal integrations can be performed by taking residues at 
branching points and at the zeros of the polynomial M{p). The examples of diagrams are collected in 
the next section. 

5.2 Low-order corrections 

We begin with presenting several low-order corrections to the free energy. 

-1- -5^0,1 = i~^ 

(87) 

This gives 

Jc^4^ y{Q) 27ri JcM 2m 

i.e., we obtain the semiclassical Dyson term |15| . 

Next two diagrams cannot be presented in the free-energy form, so we present the expressions for 
the one-loop resolvents: 

Wifi{p)= p-^O 
We also present two first "regular" terms of the free-energy expansion: 

.^ (90) 
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2 • ^2,0 = 2 <VkO +2 <{ ) + O-^^'^O 

V^ (91) 

The free-energy diagrammatic terms corresponding to resolvents H88|) and H89|) vanish (for ()88() see 
[Zj). The free-energy term J^i^o is the subleading correction in IMM calculated in jOj. It therefore 
remains only to calculate ()89|) subsequently integrating it to obtain the corresponding free-energy 
contribution ^0,2 • 

6 Calculating JFq 2 

We begin with demonstrating that the diagrammatical expression for ^0,2 of the form 

2 %}^^J + 2 <(j(j + V ';f-M^w J 

vanishes. For this, we observe that the first two terms are 

fy Jc>» y(9) Og "^^'2ny{(y 

and we can integrate by parts to set the derivative w.r.t. q on the exceptional vertex. To do this, 
however, we must first interchange the order of contour integrations (as we have the integral of y inside 
the outer integration in q). This interchanging yields the additional term 

and, upon integrating by parts, we have 

dqy{q) y'^iq) 

The constant part does not contribute, while for the second part we introduce the notation 

y (Q) -' 

The remaining contribution can be then graphically depicted as 

2 ^/>-^'~* 

' - (93) 

where we assume now that the order of appearance of vertices (from left to right) implies the contour 
ordering for the corresponding integrations. 
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As concerning the third term, if we collapse the outer integration (the one with the integral in y 
term) to the support P, it gives zero upon integration because it contains no singularities, and we 
remain with two terms appearing when passing through two other integration contours. Both these 
terms are of the same form as ()93p except they both come with the factors —1 and they differ by the 
integration order. Therefore, combining with (|93|1 . we have the contribution 

which can be evaluated by taking the residue aX q = £, (doubled because of two residues at two sheets) . 
It gives the integral over 2? with the integrand 2 J y ■ (y')^ • y^"^ Together with the first residue term 
(jUJ), it gives 

2 ^ dg 4 2/ fy'iq) y'{q) d fy'iq)\\ _ ^ / dq ^g ^d ( y'{q) 



■v 



d (y'iQ)\\ _ o / dq_ g d_(; 
dq\y{q))) /c(?)27r^^?^ dq\i 



(^^ 2iTiy{q) \y{q) y{q) dq \y{q) J J Jcf^) Ini''^ dq \y'^{q) 



d.q y'jq) 



and, integrating by parts, we just obtain 2 £,(,) f^^-pr = 2n, i.e., the constant independent on the 
potential. 

We can now make a guess for the actual J^o^2- It seems very plausible to expect it to be of the 
Polyakov's anomaly form ff R-^R, where R is the curvature and 1/A is the Green's function for the 
Laplace operator, which in our case is the logarithm of the Prime form. The curvature is expressed 
through the function y as R ^ u' /u- That is, we have two natural candidates for J-'o,2' 

:Fo,- [[dqdpy^lo,Eip,q/^^^ 



y{q) ' y{p) ' 

where E is the prime form, or 



^c 



0,2 



/ / dqdp\ogy{q)B{p,q)\ogy{p), 



but neither of these expressions is well defined. The first one develops the logarithmic cut at p = q and 
cannot be written in a contour-independent way; moreover, both these expressions are divergent when 
integrating along the support. We therefore must find another representation imitating this term. A 
good choice is when we integrate by part only once, we then have the expression of the form 

f dq^y_^ f dEg,g{p) log y{p)dp, (94) 

/c(9) 2711 y{q) Jo 

where the second integral is taken along just the eigenvalue support D. 

Variation of the logarithmic term in ()94|) can be presented already in the form of the contour 
integral. In what follows, we also systematically use that 

d(l ^^ ,^y'{(i) 



/ dqB{p,q)\ogy{q) = - <f , ^ -^dEq^q{p) 



yio) 



for the point p outside the contour C^ . 
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Action of Qy,j.\ on "outer" combination y'{q)/y{q) gives dqB{r,q)/y{q), which, upon integration 



dV{r) 

by parts, gives 



^ 1 / j^ /„^y'(p) 



= / dEg^gir)^ i dEp^piq)^ - i dEg^g{r)-^^ i dEp^p{q), (95) 



-"D ^ \ i / <^jy ^ \i / l^^ 



and in the second term we recognize the two first diagrams of (|89)) (with factors —1). 

Action of Qy/^\ on "inner" logy(p) gives B{r,p)/y{p), and in order to correspond to our diagram- 
matic technique we must interchange the order of integration w.r.t. q and p. And for the last action of 
Qyf^\ on dEp^p{q), we use (|6ip . That is, implying the contour ordering as above (from left to right) and 
indicating explicitly all the y-factors in the vertices, we graphically obtain for the last two variations: 




y 



+ 



1 

y 

p 



logy 



p 




IL 
y 



v. 
<L y 

y' 1/ y' 

^IF -\- r y^ + r ^ 

s p s 



y 



-2 r 1 ^^ - r .- ^ 

s py 



-9j^ — ty-\ - r. ^dJX 



p P\ yip) y I s ^ \v 



s 7?I 




^. 1 5 1, y. 

p y{p) P\ y 



r 



i. — ^yi _ , — ^rc-< 



py s\y^ 
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The first line, together with the second term in H95|) . just give (with the minus sign) the desired term 
VFo.2(?") H89() . the first term in the second line is canceled by the first term in ()95() . so the only mismatch 
is due to the second term in the second line. This term acquires the form 

cW y{s)\y{s) y^{s) J \ds \y{s) 

where we use the standard notation of [Jj (cf. (|1H|) ). But this term can be already calculated just by 
taking residues at the branching points, which gives (for the comprehensive calculation, see jlUj ) 

E 3?... (2,737^) = -,wv) '°« (n *«'A(.) j , (%) 

where Ma = M{p)\ are the first moments of the potential and A(;u) is the Vandermonde deter- 

minant in ^Q,. 

Therefore, combining all the terms, we conclude that 

^"'- = - 4. ^'mlo •"''•'^^ "^ '<"'* - k ""- (n "^''^f")) . 07) 

and we see that, in analogy with the answer obtained by Wigmann and Zabrodin j3Uj for the corre- 
sponding correction in the normal matrix model case, we have quantum correction term similar to the 
one in Tifi (the second term in (jHZI))- 

This completes the calculation of the exceptional term .7^3,2 and we therefore have all terms of the 
asymptotic expansion of the eigenvalue model (0) . 
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